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ABSTRACT 

The evolution of protostellar outflow is investigated with resistive magneto- 
hydrodynamic nested-grid simulations that cover a wide range of spatial scales 
(~ lAU-lpc). We follow cloud evolution from the pre-stellar core stage until the 
infalling envelope dissipates long after the protostar formation. We also calculate pro- 
tostellar evolution to derive protostellar luminosity with time-dependent mass accre- 
tion through a circumstellar disk. The protostellar outflow is driven by the first core 
prior to protostar formation and is directly driven by the circumstellar disk after pro- 
tostar formation. The opening angle of the outflow is large in the Class stage. A large 
fraction of the cloud mass is ejected in this stage, which reduces the star formation 
efficiency to ~ 50 %. After the outflow breaks out from the natal cloud, the outflow 
collimation is gradually improved in the Class I stage. The head of the outflow travels 
more than ~ 10 5 AU in ~ 10 5 yr. The outflow momentum, energy and mass derived 
in our calculations agree well with observations. In addition, our simulations show the 
same correlations among outflow momentum flux, protostellar luminosity and enve- 
lope mass as those in observations. These correlations differ between Class and I 
stages, which is explained by different evolutionary stages of the outflow; in the Class 
stage, the outflow is powered by the accreting mass and acquires its momentum from 
the infalling envelope; in the Class I stage, the outflow enters the momentum-driven 
snow-plough phase. Our results suggest that protostellar outflow should determine the 
final stellar mass and significantly affect the early evolution of low-mass protostars. 
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1 INTRODUCTION 

Molecular outflows are ubiquitously observed in the star 
forming region, which indicates that young protostars gen- 
erally drive the outflows. The molecular outflow can dump a 
large fraction of cloud matter into the interstellar space, and 
only the remaining gas around the protostar contributes to 
protostellar mass growth. Therefore, the molecular outflow 
controls the resulting stellar mass and significantly affects 
the star formation process. The star forms in a gravitation- 
ally contracting cloud. Although the specific outflow driving 
mechanism is uncertain, the molecular outflow, in principal, 
is powered by the gravitational energy of the infalling mat- 
ter released in the gravitationally contracting cloud. The in- 
falling matter, or infalling envelope, exists onl y in the early 
phase of star formation (Class and I stages; lAndre et al.l 
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1 19951 : [Andre fc Montmerldl 19941 ) during which powerful out- 
flows are often observed. Thus, observation of the molecular 
outflows provides a clue for understanding the early phase 
of star formation. 

Since the f irst discovery of molecular outflow 
|Snell et all Il980h . more than 300 outflows have been 
observed in various s tar forming region s dWu et al. 112004 
lHatchell et "all 120071 '). ICabrit fc BertoutJ (| 19921 ') found that 
the outflow momentum flux with 16 outflow samples sys- 
tematically increases w ith the stellar bolometric luminosity. 
iBontemps et al.l (| 19961 ) showed that the outflow momentum 
flux also correlates well with the (infalling) envelope mass 
for Class and I stages. Moreover, they also argued that 
the outflow power decreases with time during the accretion 
phase and that the outflow properties qualitatively differ 
between Class a nd I protostars . Howe ver, with consider- 
able data scatter, lHatchell et al.l (|2007h failed to confirm 
that Class I protostars generally have a lower momentum 
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flux than Class sources. Recently, ICurtis et al.l (|2010l ) 
analysed the outflow properties of 45 samples and reported 
a decrease in outfl ow momentum flux fro m the Class to I 
stage, as shown in iBontemps et al.l (|l996l ). 

Powerful outflows are f requently observed near the 
youn gest (Class 0) objects |Bachiller fc Gomez- Gonzalez! 
1 19921 ), indicating that vigorous outflow emerges in the very 
early evolutionary phase in which th e protostar has attaine d 
a small fraction of its final mass (jBontemps et all 1 19961) . 
These observations suggest that molecular outflow affects 
the early evolution of newly born stars. However, we cannot 
directly observe the outflow driving region, which is deeply 
embedded in a dense infalling envelope. Thus, it is difficult 
to specify the outflow driving mechanism with only obser- 
vational results. Theoretical modelling or numerical simula- 
tions are necessary to understand and thereby resolve this 
issue. 

Sinc e the discovery of well -collimated jet-like flows (op- 
tical jets: iMundt fc Friedlll983l ), it has been postulated that 
low- velocity wide-angle flows, or molecular outflows, are en- 
trained by these jets. This simple notion has been preva- 
lent because both high-velocity jets and low-velocity out- 
flows are comprehensively explained. Many authors have 
proposed various entrainment models llCabrit et ah 1 1 19971 ; 
iRicher et aLll2000l ; iLee et al-llioool ; lArce et al. 1120071 ) to an- 
alytically or numerically study the outflow driving mecha- 
nism in which the jet is artificially injected into the ambient 
medium to entrain the infalling material. However, it is dif- 
ficult to specify the outflow driving mechanism because an 
abundance of free or unknown parameters is available to re- 
produce low- velocity outflow entrained by high-velocity jets 

(953- 

Conversely, a completely different the oretical concept o f 
the molecular outflow has been proposed. iTomisakal (|2002l ) 
calculated the gravitational collapse of a molecular cloud 
core and demonstrated that both high-velocity and low- 
velocity flows are launched from different obj ects; the low- 
velocity flow is launched fro m the first core (|Larson|[l969l ; 
iMasunaga fc Inutsukal I2OO0I ) . and the high- velocity flow is 
launched near the protostar. This concept has been sup- 
po rted by many oth er cloud-collapse simulations (see review 
of iMachiddlirjlldh . Therefore, in such works, the authors 
followed the evolution with the absence of free parameters 
to control the jet and outflow. However, they could not cal- 
culate long-term ev olution of outflow because the protosta r 
itself was resolved (|Tomisakal |2002| ; iMachida et al 1 l2008bh . 
and the numerical timestep be came increasingly shor t as the 
protostar and jet evolved (e.g. IMachida et aL 12008a '). 

Long-term cloud co llapse simulations were condu cted 
in our previous work (|Machida fc Matsumotd |2012| ) us- 
ing the sink-cell technique. We demonstrated that the 
first core evolv es to the c i rcums t ellar disk afte r proto - 
star formation (|Batd 1 19981 . |2010| ; IMachida etaLI l2010al ). 
The low-velocity flow, which is launched from the first 
core prior to the protostar formation, is driven by the 
circumstellar disk after its formation. Very recently, out- 
flow driven by the first core (candida t e) were observed 



by several authors jDunham et al.l 2011 ; Chen et al. Il2010l ; 
lEnoch et all |2010| ; iPineda et all 1201 ll ; lchen et al l \20m . 

A considerably young outflow was also observed around 
an extremely young prestellar core jTakahashi fc Holl2012l ; 



iTakahashi et al.l 120121 ) . These observations agree well with 
our simulation results. 

In this paper, we extend our previous work to directly 
compare outflows simulated in the later evolutionary stages 
with observations by calculating cloud evolution from the 
prestellar cloud stage until the protostar evolves into the 
Class I or II stage. We compare the resulting outflow prop- 
erties such as the momentum flux, energy and shape with 
observation data. In addition, the protostellar evolution is 
numerically calculated with time-dependent accretion his- 
tories obtained in the simulations. This process enables 
the examination of the correlations among outflow momen- 
tum flux, protostellar luminosity an d envelope mass, which 
are suggested by obser vations (e.g. IBontemps et al.l 1 19961 : 
ICabrit fc Bertoutlll992l ). 

This paper is structured in the following manner. The 
framework of our models and the numerical method are de- 
scribed in §2, and the numerical results are presented in §3. 
We compare calculation results with observations in §4, dis- 
cuss the model parameters and outflow models in §5 and 
summarize our results in 86. 



2 NUMERICAL METHOD 

In this study, we calculate the evolution of both cloud cores 
and protostars. With our magnetohydrodynamic (MHD) 
code, we calculate the evolution of a collapsing core from 
the prestellar stage until the point at which the collapsing 
core, or the infalling envelope, dissipates after the protostar 
formation. In our simulations, we demonstrate that proto- 
stellar outflows are driven by circumstellar disks. To fol- 
low the long-term evolution over the entire main accretion 
phase, we adopt sink cells for masking the specific vicinity 
of the protostar. In addition, we calculate the protostellar 
evolution with time-dependent accretion histories obtained 
in the MHD simulations. This combination of MHD simula- 
tions and stellar evolution calculations enable us to link the 
evolution of protostellar outflow with protostellar evolution. 
In this section, we first describe the method and settings of 
the MHD calculations and explain our modelling process of 
protostellar evolution. 



2.1 MHD Calculation 

2.1.1 Initial Settings 

To investigate the evolution of cloud cores, we use three- 
dimensional resistive MHD equation s including self-gravity . 
The numerical method is described in Machida ct al. 
IMachida etHI (|2005al ) and IMachida et al. I (|2005b1 ). The 
basic equations, resistivity and sink cell treatment are th e 
same as those reported in IMachida fc Matsumotd (|2012l ). 
In the calculation, instead of solving the energy equa- 
tion, we use a barotropic e quation of state (eq. [5] of 
IMachida fc Matsumotol |2012| ) which mimics the thermal 
evolution of the collapsing prest ellar cloud core (|Larsonl 
1 19691 ; IMasunaga fc Inutsukal |2000| ). This treatment makes 
a long-term calculation possible. It should be noted that, 
however, the barotropic equation of state is not ade- 
quate especially long after the protostar formation because 
the protostellar luminosity can heat up the surrounding 
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gas dMasunaga fe Inutsukall200d ; IWhitehouse fc Batdl2006l ; 
lKrumholj|2006l ). 

As the initial state, we assume an isolated cloud core 
embedded in an interstellar medium. We adopt a spheri- 
cal cloud with a critical Bonnor-Ebert (BE) density pro- 
file, pbe, in which a uniform density is adopted outside the 
sphere (r > R c i, where R c i is the critical BE radius) to mimic 
the interstellar medium. We prohibit gas inflow at r — R c to 
strictly avoid mass inflow from outside the core, while we do 
not prohibit mass outflow at the boundary between the BE 
sphere and interstellar medium. Thus, the outflowing gas 
can freely escape from the BE sphere through protostellar 
outflow. Hereafter, we refer to the gravitationally bound gas 
cloud within r < R c as the host cloud. It should be noted 
that we confirmed that the total mass of the host cloud is 
well conserved during the calculation before the protostellar 
outflow reaches the cloud boundary. 

Because the critical BE sphere is in equilibrium, we 
increase the density by a factor of / = 1.68 to pro- 
mote contraction, where / is the density enhancement 
factor that represents the stability of the initial cloud. 
The cloud stability is generally represented by a param- 
eter «o (= Et/Eg), which is the ratio of thermal en- 
ergy (E t ) to gravitational energy (-Eg). The density en- 
hancement factor of / = 1.68 corresponds to qq = 0.5 
(|Machida et all 120061 ; iMachida fc Matsumotd 12012] ) . The 
density p rofile of the initial cloud is th e same as that re- 
ported in IMachida fe Matsumotd |20I2i ). 

For a dimensional BE density profile, we adopt an 
isothermal temperature of T = 10 K and a central num- 
ber density of n c — 6 x I0 5 cm~ 3 . With these parameters, 
the critical BE radius (tbe) is tbe = 6.1 x I0 3 AU. For 
typical models, the cloud core has a critical BE radius of 
Rd = tbe, where R c i is the cloud core radius. The mass 
inside r < R c \ for the typical models is M a \ = 1.05 Mq. 
In addition, to investigate the mass dependence of the host 
cloud on the evolution of protostellar outflow, we create two 
exceptional models with different initial cloud mass, which 
changes the cloud radius R c . These models are 1.5 and 
2 times the critical BE radius; the cloud radius for these 
models is R c = 1.5 tbe (= 9.2 x 10 3 AU) and Rc = 2r BE 
(=1.2 x 10 4 AU), respectively. The initial cloud mass for 
these models is M c \ = 1.6 Mq for the model with r = 1.5 R c \ 
and Mci = 2.1 Mq for that with r = 2 Rd, respectively. The 
model name, host cloud radius and mass are listed in Ta- 
ble ED 

In each model, the cloud rotates rigidly around the z- 
axis in the r < R c \ region and a uniform magnetic field 
parallel to the z-axis, or rotation axis, is adopted in the 
entire computational domain. We parameterized the initial 
magnetic field strength and rotation rate. The magnetic field 
strength is scaled using the central density po and sound 
speed c s ,o as 

& = B 2 /(47rpo4o), (1) 
and the rotation rate is scaled using the central density as 
a; = fi /(47rGpo) 1/2 . (2) 

With these parameters, we created ten models as listed in 
Table [T] The dimensional magnetic field strength (Bo) and 
angular velocity (Qo) for each model are also described in 
Table [T] To simply characterize models, we estimated the 



ratios of rotational and magnetic energies to the gravita- 
tional energy, /3 (= E rot /E gISLV ) and 70 (= Ema g /E grav ), 
where E rot and £7 ma g are rotational and magnetic ener- 
gies, respectively, and summarized them in Table [T] We also 
estimated the mass-to-flux ra tio M/<E> of each host cloud. 
iMouschovias fc Spitzerl (|l976l ) derived the following critical 
mass-to-flux ratio 

where the constant £ = 0.48 (jTomisaka et al. Ill988al lbh. The 
mass-to-flux ratio normalized by the critical value p is de- 
scribed as 

Models include a normalized mass-to-flux ratio of 1.8 ^ p ^ 
23. The normalized mass-to-flux ratio is also listed in Ta- 
bled] 



2.1.2 Sink Cell and Numerical Method 

To realize the long-term calculation of star formation, we 
adopt a sink at the centre of the cloud. The detailed proce- 
dure fo r intr oducing sink cells is described in IMachida et all 
(|2010ah and lMachida fc Matsumotd i|2012f ). Here, we briefly 
describe the process. We begin the calculation without a sink 
and calculate the cloud evolution for the pre-stellar gas col- 
lapse phase without a sink. Then, we identify the protostar 
formation in the collapsing cloud when the number density 
exceeds n > nthr = 10 13 cm~ 3 at the cloud centre. After 
protostar formation, in the region r < r s i n k = 2AU, the 
gas with a number density of n > 7i t hr is removed from the 
computational domain and added to the protostar as a grav- 
itating mass in each time step. Thus, for each time step, the 
accretion mass onto the protostar is calculated as 

M acc = / [p(i, j,k) - pthv]dV. (5) 

We store the mass accretion rate in one-year increment. It 
should be noted that because simulation timestep is shorter 
than At < 0.01 yr, the mass accretion rate used in this study 
is averaged over > 100 timesteps. 

To calculate the outflow driving region and long- 
distance propagation of outflow, a wide range of spatial scale 
from ~ 1 AU (the scale of the circumstellar disk) to ~ 1 pc 
(the scale of the evolved protostellar outflow) must be cov- 
ered. To resolve such considerably different scales, we use the 
nested grid method (|Machida et al. |[2005al lbh. Each level of 
a rectangular grid has the same number of cells (64 x 64 x 32) 
and the grid size and cell width is halved for every increase 
in grid level. The calculation begins with five grid levels 
(I — 1 — 5) . The fifth level of the grid (I = 5) has a box size of 
is = 2 tbe = 1.2 x 10 4 AU. Thus, the host cloud for models 
1-8 is just embedded in the I = 5 grid, while that for models 8 
and 9 is embedded in the / = 4 grid. The first level of the grid 
has a box size of Li = 2 5 r BE = 2.0 x 10 5 AU (=0.97 pc) filled 
with low-density interstellar medium pism = O.Olpo outside 
r > R c i, where pism is the density of the interstellar medium. 
Thus, we are able to calculate the propagation of the proto- 
stellar outflow in the region of < 2.0 x 10 5 AU. However, the 
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protostellar outflow never reached the computational bound- 
ary by the end of the calculation in any of the models. After 
the calculation begins, a new fi ner grid is generated b efore 
the Jeans condition is violated |Truelove et al. 1 1 19971 ') . We 
set the maximum grid level of I = 12 with a box size of 94 AU 
and cell width of 1.46 AU. Thus, we are able to resolve the 
structure from ~ 1 AU to ~ 1 pc. We calculated the cloud 
evolution with various spatial resolutions, including various 
cell widths and grid sizes, to investigate the required spatial 
resolution for outflow driving. With these calculations, we 
checked the convergence of outflow momentum and energy 
and confirmed that the spatial resolution adopted in this 
study is sufficiently high to investigat e the evolution of the 
outflow (|Machida fc Matsumotoll2012h . 



where the first and second terms on the right-hand side rep- 
resent the stellar and accretion luminosities, respectively. 
We refer to the total luminosity L tot simply as protostellar 
luminosity hereafter. In all of the examined cases, the ac- 
cretion luminosity I/ acc dominates the stellar luminosity L* . 
Because the protostellar evolution advances in the accretion 
timescale M/M in this case, we update the stellar mod- 
els each time the stellar mass increases by ~ 1 % with the 
accretion rate averaged over the timestep. The very short 
accretion variability is eliminated through this procedure, 
however, the stellar structure change are minimal over such 
short periods. 



2.2 Protostellar Calculation 

We also calculate protostellar evolution with variable mass 
accretion histories obtained through the MHD simulations. 
We numerically solve the stellar structure equ ations by con- 
sider i ng the effects of the m ass accretion (e.g. IStahler et alj 
ll98d ; |Palla fc Stahlej[l99lh . Our numerical codes have been 
developed in our previous work to examine high-mass and 
low-mass protostellar evolution with various accretion histo- 
ries ( e.g. iHosokawa fc Omukail [20091 ; iHosokawa et alj|20ld 
l201lh . In this paper, we calculate the protostellar evolution 
separately from the MHD simulations, which give us ac- 
cretion histories resulting from the interplay between the 
powerful outflow and infalling envelope. 

Our evolutionary calculation begins with a tiny initial 
model with a mass of 3 x 10~ Mq . We follow the evolution 
in which the stellar mass increases with a provided accretion 
history. Mass accretion could significantly affect the stellar 
interior structure. An important and unknown quantity of 
the mass accretion is the thermal efficiency, i.e. the specific 
entropy of gas settling onto the stellar surface. In the stel- 
lar evolution calculations, the accretion thermal efficiency is 
controlled by the outer boundary conditions of the models. 
In this paper we adopt the shock outer boundary condition, 
which supposes that the spherical accretion flow directly hits 
the stellar surface to form an accretion shock front (e.g. 
IStahler et al.|[l980l ). This condition corresponds to a rela- 
tively hot mass accretion because part of the entropy gener- 
ated at the accretion shock front is efficiently absorbed into 
the stellar interior. Some authors argue that, in the low-mass 
star formation with an accretion rate of < 10 -5 Moyr -1 , 
the accretion should be significantly colder. However, the ex- 
act value of the accretion thermal efficiency remains highly 
uncertain, and this issue is beyond the scope of our study. 
With the shock boundary, protostellar evolution is nearly 
independent of arbitrary initial models with unknown prop - 
erties (such as mass and radius; e.g. IStahler et al.l fl980h . 
With the outer boundary of the cold accretion, however, 
protostellar evolution differs wit h various initial mode ls as 
well as accretion histories (e.g. iHartmann et al.lll997T ). In 
this paper, we adopt the shock boundary to focus only on 
variations of protostellar evolution with various accretion 
histories. If the mass accretion was colder than the assumed 
value, the resulting stellar radius would be relatively smaller 
than that determined by our results. 

The protostellar calculations provides the evolution of 
the total stellar luminosity 



3 RESULTS 

3.1 Classification of Evolutionary Stage 

We calculated the cloud evolution for ten models with dif- 
ferent cloud parameters, as listed in Table [T] To characterize 
the evolution and to compare our results with observations, 
we defined three different evolutionary stages, Class 0, I, and 
II stages, using the envelope mass as 

Class stage : (M env /M cnv ,o) ^ 0.5, 
Class I stage : 0.1 < (Af cnv /M onv , ) < 0.5, 
Class II stage : (M onv /M env , ) < 0.1, 

where M e nv,o is the initial cloud mass, and the envelope mass 
A/ onv is defined as 

M anv = Mtotai (r < Rci) - Afdisk, (7) 

where Mtotai (r < Rci) is the total mass in the region of 
r < Rd , and Mdisk is the mass of the rotation -supported disk 
identifi ed accor ding to the prescription in | Machida et al.l 
(|2010ah (see also lMachida fc MatsumotolbOli T Prior to the 
formation of the protostar and rotating disk, the envelope 
mass Menv coincides with the initial cloud mass M en v,o- It 
should be noted that the protostellar mass M ps is not in- 
cluded in equation j7} because the gas composing the pro- 
tostar is removed from the computational domain, as de- 
scribed in t)2.1.2l In addition, the mass of the host cloud 
with a radius of r = -Rci (see t)2.1.1|) differs from the enve- 
lope mass by the mass of the rotation-supported disk. The 
gas falling onto the sink converts into a protostar and is re- 
moved from the computational domain, and part of the gas 
is expelled from the host cloud by the protostellar outflow. 
Therefore, the envelope mass never increases but decreases 
with time. 

In this paper, we defined Class 0, I and II stages by 
using only the env elope mass. With observational results, 
lAndre et al] (|l993l ) originally defined Class protostars as 
objects which have M ps /M ollv < 1. This indicates that the 
envelope still retains approximately a half of the initial cloud 
mass during the Class stage. Thus, we defined the tran- 
sition period between the Class and I stages at the point 
at which 50% of the initial cloud is accreted onto the pro- 
tostellar system (protostar plus the circumstellar disk). It 
should be noted that, in our definition, the envelope mass 
includes the outflowing gas in addition to the infalling gas. 
We can determine this transition period using only the mass 
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of the infalling gas. However, we believe that our definition 
is plausible for comparing calculation results with observa- 
tions because it is difficult to separate the infalling gas from 
the outflowing gas in observations. 

The use of only envelope mass creates difficulties in 
determining the transition period between Class I and 
II stage because these classes are observationally deter- 
mined through spectral energy distribution. However, such a 
method is useful in defining the classes with a single param- 
eter (M env ). In this paper, we focus mainly on the evolution 
during the Class and I stages, and briefly comment on 
the Class II stage. In addition, it is considered that nearly 
all of the infalling envelope has already disappeared in the 
Class II stage. Thus, we roughly define the transition pe- 
riod between the Class I and II stage as the period at which 
the envelope mass reaches 10 % of the initial cloud ma ss, as 
previously described (see also I Vorobvov fc~Ba su 2006). 

Table [2] summarizes the calculation results for each 
model. The mass of the protostar (M ps> o), rotating disk 
(-Mdisk,o) and protostellar outflow (M ou t,o) at the end of the 
Class stage are listed from the second to fourth columns, 
respectively, of Table [2] The protostellar luminosity (£ p8 ,o) 
at the end of the Class stage is listed in the fifth column. 
The duration (to) of the Class stage after protostar forma- 
tion is described in the sixth column. The same quantities 
(M ps ,i, Mdisk.i, M ut, i, i ps ,i and ti) at the end of the Class 
I stage are listed in the seventh to eleventh columns. These 
quantities are not described for models 1, 4, 5, 9 and 10 be- 
cause the calculation did not proceed in the Class II stage 
in such models. The mass of the protostar (M p8j n), rotating 
disk (MdiBk.n), protostellar outflow (M ou t,n) and envelope 
(Menv) at the end of the calculation are listed in the twelfth 
to fifteenth columns. In this section, we describe the evolu- 
tion of the cloud and protostellar outflow for a typical model 
(model 3) in i]3.2l and we summarize the results for all mod- 
els in JO] 

3.2 Evolution of Protostellar Outflow for Typical 
Model 

3.2.1 Structure and Collimation of Protostellar Outflow 

In this subsection, we show the cloud evolution for a typical 
model (model 3). At the initial state, the host cloud has a 
radius of r = 6.1 x 10 3 AU and a mass of M = 1.05 Mq. 
The initial magnetic field strength and rotation rate are 
Bo = 25^G and fi = 1.0 x 10~ 13 s _1 , respectively. The 
mass-to-flux ratio is /J, = 7. Observations show that molec- 
ular cloud cores have a mass-to- flux ratio of 0.8 <; <^ 7.2 
with a median value of A ~ 2 (e.g. ICrutcherlll999l ). Thus, 
Model 3 has a somewhat weaker magnetic field than that 
of the observed values. The initial rotational energy of the 
cloud is ft = 0.01. Observations show that molecular cloud 
cores have 10~ 4 < /3p < 1.4 with a ty pical value of /3o ~ 0.02 
(|Goodman et al. Ill993l ; ICasellilbOOa . e.g.). The adopted ro- 
tational energy is also slightly smaller than that of the ob- 
servational estimates. 

Figure [1] depicts several key objects for guiding the evo- 
lution of protostellar outflow in various spatial scales. In 
this paper, we used three different elapsed times: the elapsed 
time after the cloud begins to collapse (t), that after proto- 
star formation (t ps ) and that after outflow emergence (tout). 



Figure [TJ shows the structure at t = 2.178 x 10 5 yr (Class 
II stage), which corresponds to t pB — 1.447 x 10 5 yr and 
tout = 1.473 x 10 5 yr. The protostar formation epoch t pSj o 
and outflow emergence time t ou t,o are listed in Table [3] 

As indicated in Figure [TJi, the protostellar outflow has a 
size of ~ 2.5 x 10 4 AU on one side at this epoch. Because the 
radius of the host cloud is R c \ = 6.1 x 10 AU, the vertical 
length of the protostellar outflow is approximately four times 
the host cloud radius. However, the horizontal width of the 
protostell ar outflow is comparable to th e host cloud radius 
(Fig. ITb). iMachida fc Matsumotol (|2012h reported that the 
width of the outflow reflects the host cloud radius because 
the protostellar outflow is anchored by the cloud-scale mag- 
netic field lines and can widen up to the cloud scale. Thus, 
we can observationally determine the size and mass of the 
core from the (maximum) width of the protostellar outflow. 
The blue sphere in Figures [TJi and b corresponds to the host 
cloud. In Figures [if) and c, inside the host cloud, the dense 
infalling envelope has a torus- like structure and can be iden- 
tified by the orange iso-density surface of n — 7 x 10 4 cm -3 . 
A pseudo-disk that is not rotation-supported disk is evident 
within the dense infalling envelope and is represented by 
the green iso-density surface. As indicated in Figures [TJc, d, 
and e, the protostellar outflow is launched from a rotating 
disk that is surrounded by the pseudo-disk in the very cen- 
tre of the core. At this epoch, the rotating disk has a size 
of ~ 100 AU, and the mass of the protostar and rotating 
disk are M ps = 0.45 Mq and Mdisk = 0.11 Mq, respectively. 
As shown in Figure [ljs, the rotating disk is vertically pene- 
trated by the magnetic field lines that are significantly in- 
clined against the disk rotation axis corresponding to the 
z-axis. Figure [TJi shows that inside the protostellar outflow, 
the magnetic field lines are strongly twisted due to the disk 
rotation. In summary, at this epoch, the outflow driving re- 
gion is embedded in r <J 100 AU, while the outflow extends 
up to ~ 2.5 x 10 4 AU. 

Figure [2] shows the time evolution of the protostellar 
outflow for model 3 in the host cloud scale (~ 1.2 x 10 4 AU). 
To compare the protostellar outflow with the host cloud, we 
only show the region of 0° < 9 < 90° and 0° < 4> < 90° 
inside the host cloud, where 9 and <f> are the zenith and 
azimuthal angles, respectively. Figure [2]a. shows the initial 
state for this model. The protostellar outflow appears in 
the collapsing cloud 7.047 x 10 4 yr after the cloud begins 
to collapse. The outflow evolves and maintains a prolate 
shape, as shown in Figures [2f> and c. The outflow reaches the 
boundary between the host cloud and interstellar medium 
at t ps = 2.521 x 10 4 yr after protostar formation. Then, the 
protostellar outflow penetrates the host cloud and flows into 
the interstellar space. The gas is ejected from the host cloud 
by the protostellar outflow for t pa > 2.521 x 10 4 yr. 

As shown in Figures (2js and /, the volume occupied by 
the protostellar outflow inside the host cloud gradually in- 
creases with time. The protostellar outflow propagates along 
the magnetic field line. The magnetic field lines open up as 
the distance from the equatorial plane increases and devel- 
ops an hourglass- like configuration. Thus, the opening an- 
gle of the protostellar outflow also gradually increases with 
inc reasing distance from the equ atorial plane (for details, 
see IMachida fc Matsumotol 120121 ). As a result, the proto- 
stellar outflow can sweep up and incorporate a large frac- 
tion of the infalling gas and eject it into the interstellar 
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space. Therefore, protostellar outflow reduces the star for- 
mation efficiency jNakano et al~ 19951 ; iMatzner fc McKeel 
l200d ; iMachida fc Matsumotol 120121 ). After a large fraction 
of gas in the infalling envelope falls onto the protostellar 
system, or is expelled into the interstellar space (i.e. in the 
Class I or II stages), the volume occupied by the protostel- 
lar outflow becomes small. The protostellar outflow shown 
in Figure [2]f is slimmer than that in Figure [2]2 at its root, 
or near the outflow driving region because the protostellar 
outflow weakens as the infalling envelope is depleted. Thus, 
during the Class I and II stages, the mass ejection rate from 
the host cloud gradually decreases. 

After the protostellar outflow breaks out of the cloud, 
it propagates into the interstellar space. Figure [3] shows the 
evolution of the protostellar outflow after the outflow ver- 
tical length exceeds the size of the host cloud; the scale of 
each panel is different to describe the entire region of the 
outflow. The protostellar outflow is loosely collimated and 
has a relatively wide opening angle just after it penetrates 
the host cloud (Fig. [3]a) . Then, the collimation of the out- 
flow is gradually improved with time (Figs. [3f*-<i) because the 
outflow extends only in the vertical direction and expands 
minimally in the horizontal direction. As shown in Figure [31 
the outflow always has a width comparable to the size of the 
host cloud. As a result, the well-evolved outflow is well col- 
limated (Figs. and d). Such a well-collimat ed outflow is 
often observed in the star f orming region (e.g. iHirano et ail 
l2006l : IVelusamv et al.ll2007r i. As shown in Figure EH, the pro- 
tostellar outflow for model 3 reaches up to ~ 10 AU at the 
end of the calculation. 

To investigate the morphological evolution of the pro- 
tostellar outflow, we present the shape of the outflowing re- 
gion at each epoch in Figure U As shown in the lower right- 
hand panel, before the outflow reaches the cloud bound- 
ary (t ps < 3 x 10 4 yr, epochs [1] - [4]), the outflow evolves 
and maintains nearly the same ratio of vertical to horizontal 
length. Thus, during this period, the outflow opening angle 
rarely changes, and the outflow collimation is poor. After 
the outflow reaches the cloud boundary (epochs [5] - [8]), 
the outflow extends only in the vertical direction as shown 
in the left panel of the figure. The opening angle of the 
protostellar outflow gradually becomes smaller with time. 
Therefore, the outflow has a relatively wide opening angle 
in the earliest evolutionary stage, while a well-collimated 
outflow is realized in the later stage. 

On the contrary, in a cloud scale, or a scale smaller than 
the cloud radius, the opening angle increases with time, as 
indicated by grey arrows over epochs [4] - [6] in the lower 
right panel. Because the observations usually focus on the 
dense outflowing gas that is located near, or inside the cloud 
core, the opening angle of the outflow may be obse r ved to 
apparently increase with time. IVelusamv fc Langerl (|l998l ) 
observed outflow-infall interaction in IRS1 in B5 with 12 CO 
and 18 CO emissions and showed that th e opening angle 
of the outflow should increase with time. lArce fc Sargenj 
(|2006h also reported that the outflow cavity widens as the 
envelope mass decreases with different observations. These 
observations are consistent with our results in the cloud 
scale. Because the outflow has a wide opening angle in the 
cloud scale, it effectively limits the gas accretion onto the 
protostellar system. 



3.2.2 Envelope Mass and Outflow Momentum 

The mass of the infalling envelope decreases with time as gas 
accretes onto the star and disk. Because the outflow is pow- 
ered by the mass accretion, the protostellar outflow ceases as 
the infalling envelope gets depleted. To illustrate this pro- 
cess, we show the density distribution of the infalling en- 
velope (each upper panel) and outflow momentum (p|u ut|; 
each lower panel) on the y — cutting plane in Figures [5] 
and [6] Figure [5] shows the evolution of the entire host cloud 
with a box size of 1.2 x 10 4 AU, while Figure|S]shows the evo- 
lution of the infalling envelope around the outflow driving 
region with a box size of 1500 AU. 

As shown in the upper panels of Figure [5] the den- 
sity of the infalling envelope decreases with time. At the 
end of the calculation, the density just inside the host cloud 
p(r = 6.1 x 10 4 AU) is much less than 10% of that at the 
initial state. As shown in Figures 03 and d, a dense in- 
falling envelope with 10 4 cm -3 < n < 10 6 cm -3 remains 
for t <J 1.3 X 10 5 yr, and the protostellar outflow has a rel- 
atively large momentum. As the gas density in the infalling 
envelope decreases, the outflow momentum also decreases, 
which is evident through a comparison of the panels in Fig- 
ure [S] panels c and d with those in Figure and /. In ad- 
dition, it is evident that the outflowing region is separated 
in Figures [5je and /, which occurs because the protostellar 
outflow is intermittently driven by the circumstellar disk in 
the later evolutionary stage (t > 1.6 X 10 5 yr). It should be 
noted that the outflow is continuously driven in the early 
stage (t < 1.6 x 10 5 yr; Figs. EM). At the end of the cal- 
culation, the envelope mass becomes less than 10% of the 
initial cloud mass, and the very weak outflow is driven by 
the circumstellar disk. 

Figure [6] shows that the outflow has a considerably 
wide opening angle near the equatorial plane during the 
main accretion phase. The opening angle increases until 
t <C 1.5 x 10 5 yr and has a maximum opening angle of ~ 130° 
at this scale. Then, the opening angle shifts to decrease for 
t J> 1.5 x 10 5 yr. Because the outflow opening angle is large 
in the early evolutionary stage, the gas accretes onto the 
circumstellar disk only from the side as seen in Figure HJf) 
and c. The mass accretion rate and outflow momentum de- 
crease with the density of the infalling envelope (Fig. [6jf) . 
Then, the outflow is intermittently driven by the circumstel- 
lar disk and has a nested structure as shown in Figure [6)2 
and /. The outfl ow finally disappears as most o f the envelope 
gas is depleted (|Machida fc Matsumotoll2012l ). 

The time evolution of the envelope mass is plotted 
against the elapsed time t in Figure 08, in which the mass 
of the protostar, rotating disk and outflow are also plot- 
ted. The figure shows that the envelope mass is conserved 
for U 6.8 x 10 4 yr. At t ~ 6.8 x 10 4 yr, the first (adi- 
abatic) core, or the rotating disk, forms in the collapsing 
cloud, and the mass o f the infalling envelope begins to de- 
creas e. The first core (|Larsonlll969l ; iMasunaga fc Inutsukal 
bnnl i forms prior to the protostar for mation and is sup- 
porte d both by gas pressure and rotation (|Saigo fc Tomisakal 
l2006h . The infalling gas continues to accrete onto the first 
core until protostar formation. After protostar formation, 
the first core becomes the ci rcumstellar disk, which is mainly 
supported by the rotation dBatell 19981 . l201ll ; IMachida et al.1 
l2010al ; lTsukamoto fc Machidall201li r After protostar forma- 
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tion at tp S ,o = 7.3127 x 10 4 yr, the gas accretes onto the 
protostar through the circumstellar, or rotating, disk. Fig- 
ure 03 shows that the circumstellar disk mass dominates 
the protostellar mass for t < 9 x 10 4 yr, or approximately 
2 x 10 4 yr after the protostar formation. The protostellar 
mass then dominates the circumstellar disk, which tends to 
become gravitationally stable. 

The outflow is launched just prior to protostar forma- 
tion (tout.o = 7.047 x 10 4 yr). As shown in 33.2.11 the out- 
flow sweeps the gas in the infalling envelope. The protostel- 
lar outflow can eject a mass comparable to the protostellar 
mass into the interstellar space from the host cloud. The 
blue dotted line in Figure [TJi represents the outflowing mass 
inside the host cloud and indicates that the outflow grad- 
ually weakens for t > 10 yr. The mass of total outflowing 
gas includes the outflowing gas within the host cloud, gas 
ejected from the host cloud and interstellar gas swept by the 
outflow. The mass of the incorporated interstellar gas is less 
than 3% of the total outflowing mass because of the very 
low ambient density adopted outside of the cloud (SJ27TTTJ) . 

The masses of the protostar, disk and outflow are plot- 
ted against the envelope mass normalized by the initial cloud 
mass in Figure0, where the evolutionary stages of Class 0, I 
and II (see i|3.1[) are denoted with background colours. Fig- 
ure shows that the outflow mass inside the host cloud, 
represented by the blue broken line, shifts to decreases just 
prior to the Class I stage. In addition, the disk mass gradu- 
ally decreases during the Class I stage. On the contrary, the 
protostellar and total outflowing mass slightly increase even 
during the Class I stage. 

The mass accretion rate onto the protostar, represented 
by solid black line, is plotted against the elapsed time t in 
Figure and against the envelope mass ratio in Figure 0f, 
where the red line represents the mass accretion rate aver- 
aged over 1000 yr. Just after the protostar formation, the 
mass accretion rate is ~ 1O~ 5 M0 yr -1 and gradually de- 
creases during the Class stage. Then, early in the Class I 
stage, the averaged mass accretion rate temporally increases 
for 10 5 yr < t < 1.5 X 10 s yr, or 0.3 < Af cnv /A/ cnv . < 0.45. 
It is evident that the mass accretion onto the protostar is 
highly time variable. This time variability is due to the disk 
instability. During this period, the circumstellar disk grows 
to become gravitationally unstable and the angular momen- 
tum is transferred by the gravitational torque in the disk 
in addition to the ma gnetic effects of protostellar outflow 
and m agnetic braking (jMachida et al ] |2010al ;l Inutsuka et al.l 
The angular momentum transfer by the magnetic ef- 
fects leads to steady accretion, whil e that by disk instability 
often causes time variable accretion jVorobvov fc Basull2006l ; 
iMachida et~al]|2011ah . 

After the infalling envelope mass decreases consider- 
ably, the gas intermittently accretes onto the protostar from 
the circumstellar disk (t J> 1.5 x 10 s yr and M e nv/M env ,o <, 
0.3). Even during this phase, the gas continuously settles 
onto the circumstellar disk from the infalling envelope with 
relatively low mass accretion rates. The disk becomes grav- 
itationally unstable with this continuous mass supply from 
the infalling envelope. The gravitationally unstable disk in- 
stantaneously amplifies the gas accretion rate onto the pro- 
tosta r, which is known as the burst phase (|Vorobvov fc Basra 
120061 ). Then, the disk mass decreases, and the mass ac- 
cretion onto the protostar temporally halts after the disk 



recovers to a stable sta te, known as the quiescent phase 
|Vorobvov fc Basul I2006T ) . As shown in Figures 0- d, how- 
ever, the protostellar mass hardly increases during such 
burst-like accretion events because of the short durations. 

The evolution of the protostellar luminosity is also pre- 
sented in Figures [Tjs and d. Following protostar formation, 
the protostellar luminosity increases with the stellar mass 
and has an initial peak of L pa ~ 20Lq at M env /M env> o ~ 0.3 
in the Class I stage. The protostar luminosity then decreases 
as the accretion rate decreases. This decrease in the proto- 
stellar luminosity d uring the Class I stage may solve the 
luminosity problem jKenvon et aL 1990l; Evans et a"l]|2009l ; 
lEnoch et ail2009h . lOffner fc McKed (|201lh also investigated 
the protostellar luminosity with different main accretion 
models and pointed out that a gradual decrease of the ac- 
cretion rate can solve the luminosity problem. At the end of 
the Class I stage, the protostellar luminosity is L ~ 3 L@. 
It should be noted that the evolution of the protostellar lu- 
minosity does not show spiky features, which are seen in 
the accretion history for M env /Af cnv ,o < 0.3, because the 
short variabilities are smeared out for calculating the pro- 
tostellar evolution. However, the stellar structure would not 
be changed with resolution of this episodic accretion because 
the increase of the stellar mass during that time is minisculc. 
Figures and d also indicate that the accretion luminosity 
dominates the stellar luminosity L, for M env /M cnv .o > 0.3 
while it is comparable to L* for M env /Af env ,o <C 0.3. 

In Figures and /, the evolution of the outflow mo- 
mentum and energy are represented by the red and blue 
lines, respectively. In the figure, these quantities are sepa- 
rately presented over the entire outflowing region and only 
the inside of the cloud by solid and broken lines, respectively. 
During the early Class stage, both the outflow momentum 
and energy increase. However, just prior to the Class I stage 
the outflow momentum and energy for the inside of the cloud 
begin to decrease because the accreting matter, or the en- 
velope mass, that provides the driving force for the outflow 
is depleted during the Class I stage. At the end of the Class 
I stage, the outflow momentum and energy inside the cloud 
are approximately three orders of magnitude smaller than 
their peak values. No powerful outflow is driven by the cir- 
cumstellar disk during the Class II stage. In this stage, the 
ejected outflow just propagates to nearly retain its original 
momentum and energy. However, these quantities gradually 
decrease with time through interaction between the outflow 
and interstellar medium. 



3.3 Evolution of Protostellar Outflow in Clouds 
with Different Parameters 

3.3.1 Mass Evolution 

In this subsection, we explain cases with different cloud pa- 
rameters. Figure [8] presents the evolution of the envelope 
mass, which is normalized by the initial cloud mass, for mod- 
els 1-8 against the time after the cloud begins to collapse. All 
models show similar evolution of the envelope mass. Prior 
to protostar and circumstellar disk formation, the envelope 
mass is nearly constant. The rotating disk forms and the 
envelope mass begins to decrease ~ (7 — 8) x 10 4 yr after the 
cloud begins to collapse. The envelope mass halves, and the 
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Class stage ends approximately 10 4 yr after the formation 
of the protostar or rotating disk. 

This result indicates that the Class stage lasts for 
~ 10 4 yr. In fact, the durations of the Class stage to in 
the examined models are in the range of 2.0 x 10 4 yr < 
to < 8.9 x 10 4 yr (Table [2}. These values are comparab le to 
those observationally estim ated in lAndre et all (|l993l l and 
lAndre fc Montmerlej (|l994f ). Recent observati ons indicate 
that t he Class lifetime is longer than 10 4 yr. lEnoch et alj 
(|2009h observed many embedded protostars in Perseus, Ser- 
pens and Ophiuchus and estimated a Class lifetime of 1.7± 
0.3 x 10 5 yr with a relative n umber of Class and I sources 
(see also lEvans et all 120091 '). More recently, iMaurv et al.l 
(|201ll ) shows a Class lifetime of ~ 4 - 9 x 10 4 yr with 
observations of Aquila rift complex. Thus, the Class life- 
time derived in this study agrees well with observations. 

After protostar formation, the envelope mass decreases 



with M a 



oc t 



(M ei 



-2.5V 



as shown in Fig- 



ure [5] It is difficult to analytically derive the power of 
—2.5, or —3.5, because the envelope mass includes a non- 
negligible mass of the outflowing gas ejected from the host 
cloud (Fig. 0. However, this rapid decrease in the envelope 
mass indicates that the infalling gas dissipates in a short 
duration after the protostar formation, and the main accre- 
tion phase does not last for a lengthy period. That is, the 
duration of the Class I stage is not much longer than that 
of Class stage ([Enoch et al.ll2009l ; lEvans et all 120091 '). As 
described in Table [5J the durations of the Class I stage t\ 
are in the range of 7.1 x 10 4 yr < ti < 1.1 X 10 5 yr. 

Next, we comment on the parameter dependence of the 
evolution of the envelope mass. Figure [8] shows that the en- 
velope mass for models with slow initial rotation rates or 
weak magnetic fields (models 1,2, 6, and 7) decreases rapidly, 
while it for model with an initially rapid rotation or strong 
magnetic field (models 3, 4, 5, and 8) decreases slowly. Both 
the cloud rotation and magnetic field slow the cloud collapse 
and its e volution because they can support the cloud against 
gravity (|Scott fc Blacklll980l ; [Machida et al. Il2005ar ). Thus, 
clouds with rapid rotations or strong magnetic fields have 
relatively long lifetimes of the infalling envelope. 

However, the difference in cloud lifetime among mod- 
els is not significant. As shown in Figure [8l each cloud dis- 
sipates in approximately < 10 iff ,0 after cloud begins to 
collapse, where tg^o is the freefall timescale of the initial 
cloud. Thus, independent of magnetic field strength and ro- 
tation rate, clouds with the same mass (or same central 
density) have nearly the same lifetimes because gravity and 
pressure gradient force mainly control the cloud evolution, 
while rotation and mag netic field offer minor contributions 
([Machida et al. Il2005ab . 

The envelope gas gradually accretes onto the circum- 
stellar disk. Then, part of the circumstellar disk gas is blown 
away by the protostellar outflow, and part of it accretes fur- 
ther onto the protostar. The mass of the protostar, circum- 
stellar disk and protostellar outflow for models 2, 4, 5, 6, 
7 and 8 are plotted against the normalized envelope mass 
in Figure [5] In the figure, although models have the same 
initial cloud mass (Table [l}, the mass ratio of each object, 
which includes the protostar, disk and outflow, differs con- 
siderably. Furthermore, Figure [5] indicates that star forma- 
tion efficiency (e) strongly depends on cloud parameters b 
(magnetic field strength) and cj (rotation rate), where the 



star formation efficiency is defined as the ratio of the pro- 
tostellar mass to the host cloud mass e = M ps /M e nv,o- This 
result is expected because the outflow efficiency, which de- 
termines star formation efficiency, depends on the magnetic 
field strength and rotation rate of the cloud core. When the 
host cloud has no magnetic field, all the envelope, or host 
cloud, mass accretes onto the protostellar system without 
emergence of the protostellar outflow. When the host cloud 
has no rotation and no circumstellar disk forms, then all of 
the envelope mass falls directly onto the protostar. Thus, 
in the case with no rotation, we expect a star formation 
efficiency of e = 1. With rotation and magnetic field, the 
protostellar outflow can reduce the star formation efficiency 
down to £ ~ 0.5. 

As shown in Figure [9l for all models, the rotat- 
ing di sk a p pears prior t o pro t ostar formation dBatej 
19981. l201ll; IWalch et all l2009al ; iMachida et ail l2010al ; 



Tsukamoto fc Machida 20 111 - ). Just after protostar forma- 



tion, the rotation disk mass dominates the protostellar mass 
(|lnutsuka et aljfeoiOT ). In the figure, the disk mass dominates 
the protostellar mass until the end of the calculation only 
for model 2 (panel a), while the protostellar mass dominates 
the disk mass as the envelope mass decreases in other models 
(models 4, 5, 6, 7 and 8). The initial magnetic field for model 
2 is relatively weak (b = 0.05), and the angular momentum 
is not effectively transferred by the magnetic braking and 
outflow. Thus, for model 2, a massive disk comparable to 
the protostellar mass remains even in the Class II stage. 
Figure [5] also indicates that a sufficiently large disk already 
exists in the Class stage. Recent observation confirmed a 
large circumstellar disk in the Class stage |Tobin et al.l 
I2012T ). which is consistent with our results. 

Next, we focus on the mass of the protostellar outflow. 
In models 2 (panel a) and 8 (panel f), the mass of the pro- 
tostellar outflow is comparable to or larger than the proto- 
stellar mass during the Class 0, I and II stage. For these 
models, the initial cloud has a somewhat weak magnetic 
field and relatively rapid rotation. When the host cloud is 
strongly magnetized, no powerful outflow appears, as shown 
by model 5 (Fig. \§]p) , because the disk angular momentum 
is effectively transferred by the magnetic braking, and no 
massive disk to drive the protostellar outflow appears. In 
addition, when the host cloud has a small angular momen- 
tum, neither a sufficiently large disk nor powerful outflow 
appear (model 6, Fig. |9]i). Thus, the host cloud that has 
a moderately strong magnetic field and rapid rotation can 
drive a powerful, or massive, outflow, as seen in models 2 
and 8. However, even in models 5 and 6, the mass fraction 
of protostellar outflow is not negligible; at least > 10% of 
the host cloud mass is ejected from the host cloud by the 
protostellar outflow. 

Figure [9] indicates that, for all models except for model 
6, the outflow mass inside the host cloud has a peak dur- 
ing the Class stage and decreases during the Class I 
stage. In addition, as shown in the figure, the outflow 
appears prior to protostar formation and after disk for- 
mation for models 2, 4, 5 and 8, while it appears after 
the protostar formation for models 6 and 7. In general, 
the outflow is driven prior to protostar formation (e.g 
Tomisakal |2002|; iBaneriee fc Pudrit3 120061; [Machida et al 
2008b); iHennebelle fc Fromang||200Sl; iTomida et alJbOlOalT 



20121 ; (Puffin fc Pudritdl2009l ; iDuffin et alj|201ll : IPrice et ail 
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120121 ; ISeifried et afll2012h . The delayed emergence of the 
outflow for models 6 and 7 is due to the initial small ro- 
tation rate. It is expected that, in these models, the outflow 
appears after the disk acquires a sufficient angular momen- 
tum from the infalling envelope. 

Figure [9] shows that the total outflow mass is in the 
range of A/ out ~ 0.1 — 0.5 Mq. Thus, approximately 10 — 50% 
of the host cloud mass is ejected by the protostellar out- 
flow because the in itial cloud has a mass of M c i = 1.05 Mq. 
ICurtis et al.l (|2010t ) investigated the outflow mass in four ac- 
tive star forming regions and showed that the observed out- 
flows had masses in the range of 0.002 Mq < M out < 0.4 M Q 
with an average value of 0.09 ± 0.02 Mq for Class proto- 
stars and 0.06 ±0.03 Mq for Class I protostars. Although we 
expect that the outflow mass depends on the initial cloud 
mass ( H5.ip . our results roughly agree with the observations. 



3.3.2 Outflow Momentum 



The outflow momentum is a useful index to obser- 
vationall y identify the e voluti on ary stage of a pro- 



tostar (ICabrit & Bertout 19921; Bontenros et al.l 


199£ 


Arce & Sargent 20061; Wu et al. 12004; Hatchell et all 


2007 


Curtis et alJl2010HAndre et alj|200d). In Figure[T01 the out- 



flow momenta for models 1, 3, 4, 5, 6, 7 and 8 are plot- 
ted against the normalized envelope mass (left panels) and 
the elapsed time after the outflow appears (right panels). 
The variations with different magnetic fields and with dif- 
ferent rotation rates are separately presented in the upper 
and lower panels. The outflows momenta are in the range of 
0.02M a kms~ 1 < MV out < O.2Af kms _1 at their peak. 
ICurtis et all (|2010t ) observed 45 outflows in the Perseus 
molecular cloud and showed that the outflow momentum for 
Class I objects without a highly powerful anomaly (SVS13) 
is 0.10 ± 0.03A/q kms -1 on average. Thus, the outflow mo- 
mentum derived in our study is comparable to that of the 
observational estimates. As shown in the left panels in Fig- 
ure 1101 the outflow momentum is larger with weaker mag- 
netic fields (models 1 and 3; upper panel) or with higher 
angular momentum (models 3 and 8; lower panels). Both 
a moderately strong magnetic field and rapid rotation are 
necessary for driving a powerful outflow (see ^3.3. 1 [I . 

With the same initial rotation, the outflow momentum 
inside the host cloud begins to decrease at nearly the same 
point (M e nv/M onv ,o) ~ 0.6 even with different magnetic 
fields (broken lines in Fig. [10] upper left panel). On the con- 
trary, the duration for driving the powerful outflow depends 
on the cloud rotation rate (broken lines in Fig. 1 101 lower left 
panel). These results indicate that the total amount of the 
cloud angular momentum is strongly related to the duration 
of the powerful outflow driving. This duration controls the 
total amount of the ejected mass and outflow momentum 
(solid lines in Fig. fTO]) . 

The right-hand panels in Figure [10] show that outflow 
momentum is linearly proportional to the elapsed time t ou t 
after emergence of the outflow. This evolution is approxi- 
mately written as 

AiVout = 5 x 10~ 6 (^fj M e kms" 1 . (8) 
Inside the host cloud (r < R c i), relation ([8} can be applicable 



for tout < ift,o, which indicates that the duration of outflow 
driving is approximately equivalent to the freefall timescale 
of the host cloud. During this period, the gas is vigorously 
ejected around the circumstellar disk. After tout > iff,0j 
the outflow momentum inside the host cloud drastically de- 
creases, and the outflow intermittently appears as described 
in %HT2l 



3.3.3 Protostellar Luminosity and Mass Accretion Rate 

In addition to the outflow momentum, the protostellar lumi- 
nosity could represent different stages of protostellar evolu- 
tion. Figure [11] shows a time evolution of protostellar lumi- 
nosity and mass accretion rates onto the protostar in models 
2, 4, 5, 6, 7 and 8. The figure shows common accretion histo- 
ries among the models except for model 8; the accretion rate 
is nearly constant at M pB ~ 10~ 5 Afo yr -1 for t pB < 10 4 yr 
and gradually decreases after that time. The mass accre- 
tion becomes time-variable for t pB J> several x 10 4 yr and 
almost ceases with M pB < 10 _7 Mq yr -1 for t pB f> 10 5 yr. 
Reflecting such mass accretion histories, the protostellar lu- 
minosity gradually increases from ~ 0.1 L to ~ 10 Lq for 
t ps <C 5 x 10 4 yr. The luminosity then decreases as the mass 
accretion rate decreases for t pB J> 5 x 10 4 yr. 

Unlike models 2, 4, 5, 6 and 7, model 8 shows a highly 
time- variable mass accretion throughout the evolution. With 
the initial rapid rotation in this model, the circumstellar 
disk becomes highly gravitationally unstable, which makes 
the accretion history significantly time- variable. Figure 1121 
shows the density distribution in the disk for model 8. The 
disk shows a non-axisymmetric structure created by the 
gravitational instability, with which the angular momen- 
tum is transported. With the highly time-variable accretion, 
the protostellar luminosity is also significantly time- variable 
throughout the accretion phase (Fig. \ll\f ). If the mass ac- 
cretion appears through the highly gravitationally unstable 
accretion disk, it would be thus difficult to characterize the 
protostellar evolutionary stages simply from luminosity. It 
should be noted that the odd mode of the non-axisymmetric 
density perturbation is suppressed in the disk because we 
fixed the sink (or protostar) at t h e cen ter of the cloud as 
described in $2X21 iKratter et all (|2010h showed that non- 
axisymmetric m=l mode tends to develop and fragmen- 
tatio n occurs in the disk (see also ITsukamoto fc Machidal 
1201 ll ). In such a case, the mass accretion rate and proto- 
stellar luminosity may be somewhat different from those in 
Figure [TTV. 

Figure [13] summarizes the evolution of protostellar lu- 
minosity in the examined models against the normalized 
envelope mass (upper panel), the protostellar mass (mid- 
dle panel) and the elapsed time after protostar formation 
(lower panel). Model 8 is omitted here because its complex 
luminosity evolution (Fig. \ll\f ) obstructs easy viewing of the 
evolutionary tracks. In each model, the protostellar luminos- 
ity generally peaks in the Class stage and decreases in the 
Class I stage. At their peaks, the protostellar luminosities 
reach ~ 8 — 30 Lq. Figure [T3l also shows that stellar luminos- 
ity is higher with slower initial rotation (e.g. models 6 and 
7) and lower with weaker magnetic field (e.g. models 1 and 
2), which reflects the fact that the mass accretion rate is rel- 
atively lower with a more rapid rotation or weaker magnetic 
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field. The disk rotation suppresses the rapid gas accretion 
onto the protostar because the gas is supported by centrifu- 
gal force in the disk. On the contrary, the magnetic field 
promotes gas accretion with an angular momentum transfer 
by magnetic braking. Thus, protostars forming with slower 
rotation or stronger magnetic fields have the higher proto- 
stellar luminosity. 

The upper and lower panels in Figure [13] indicate that 
in each model, the protostellar luminosity during the (late) 
Class stage does not differ significantly from that dur- 
ing the (early) Class 1 stage. The observations also show a 
smooth transition o f stellar luminosity fro m the Class to I 
stages (e.g. Fig. 1 of lBontemps et alii 19961 ). On the contrary, 
at the same epoch, the luminosity difference is as large as 
one order of magnitude among models with different cloud 
parameters (Fig. [T3] upper and lower panels). This result in- 
dicates that, in observations, dispersion of the protostellar 
luminosity may be due to different properties of clouds such 
as magnetic field strength and rotation rate rather than to 
different evolutionary stages. 



4 OUTFLOW MOMENTUM FLUX: 

COMPARISON WITH OBSERVATIONS 

4.1 Momentum Flux vs. Bolometric Luminosity 

The correlations between the observed outflow momentum 
flux and protostellar luminosity, or envelope mass, should 
suggest an underlying relationship between stellar outflow 
activity and protostellar evolution. In this section, we ex- 
amine similar correlations from our simulations and compare 
them with the observations. The outflow momentum flux for 
each model is plotted against the bolometric luminosity in 
Figure [14] and against the envelope mass in Figure [T5l To de- 
rive the outflow momentum F, we calculated total outflow 
momentum (MV)out in the outflowing region (v r > c s ,o) 
and divided that value by the elapsed time t ont after the 
outflow appearance 

p = (MVQout = JV^out dv 
tout tout 

In the figures, the upper panels present the evolutionary 
tracks for the models, and the lower panels show images 
every thousand years. Models with the different initial mag- 
netic fields and different rotation rates are separately plotted 
in the left- and right-hand panels. The standard model 3 is 
plotted in all panels for comparison. 

In Figure 114b . the tracks first move from the left 
to the lower right, then turn around to the lower left. 
T hese evolutionary track s agree well with those expected 
in lBontemps et ail () 19961 . Fig. 5). For models plotted in Fig- 
ure !14b . the outflows have momentum fluxes (i.e. the vertical 
axes) in the range of 1CT 5 < F/(M e kms" 1 yr" 1 ) < 1CT 4 at 
the emergence time of outflow, which gradually decrease un- 
til the outflow disappears. On the contrary, the protostellar 
luminosities (i.e. the horizontal axes) increase after protostar 
formation and begin to decrease after reaching peak values 
of L ps ~ 10 Lq approximately 10 4 — 10 5 yr after protostar 
formation, as shown in Figure 1111 Reflecting both evolution 
of the outflow momentum and protostellar luminosity, the 
momentum flux evolves along the arrows in Figure [14b . 



In Figure [T4"b . the evolutionary track of model 8 qualita- 
tively differs from that shown in Figure [Tib . For this model, 
the angular momentum is mainly transferred by gravita- 
tional torque (Fig. I12[) . and a highly time- variable accretion 
is realized f £13.3.31 and Fig. Ill[) . Thus, with a time- variable 
accretion, the protostellar luminosity also shows a high time 
variability. As a result, model 8 shows a zigzag evolution- 
ary track of outflow momentum flux. In addition, compared 
with the models in Figure 114b . the evolutionary track in 
model 6 (and 7) differs relatively. The outflow appears long 
after the protostar formation in model 6 (and 7), while it 
appears before protostar formation in other models (Ta- 
ble [3]). For model 6, the outflow appears 1.316 x 10 4 yr after 
the protostar formation. By this epoch, the protostar suffi- 
ciently evolves and has a luminosity of ~ 30 Lq. Thus, in 
the diagram, the outflow momentum flux suddenly appears 
at (Lboi, F) = (~30L Q , ~6xl0" 6 M kms" 1 yr" 1 ). Then, 
the momentum flux moves only to the lower left because the 
protostellar luminosity has already passed the peak. 

Next, we focus on the parameter dependence of the 
evolutionary track. The upper panels in Figure Q3] indicate 
that the evolutionary track for models with stronger mag- 
netic field tends to lie in the lower right area. This occurs 
because models with stronger magnetic fields have smaller 
outflow momentum (Fig. 1 10 p but higher protostellar lumi- 
nosity (Fig. I13[) . as described in £13.31 In addition, because 
models with smaller rotation rates have a smaller outflow 
momentum and higher protostellar luminosity (JO}, their 
evolutional track lie in the lower right area. 

In Figure im with lower panels one can roughly compare 
observations with calculation results. Since the momentum 
fluxes every thousand year are plotted in the panels, it is 
expected that protostars are frequently observed in areas 
with densely grouped circles, while they are rarely observed 
in sparsely-grouped a reas. With observed Class I sources, 
iBontemps et al.l (|l996h derived the best fit for linear correla- 
tion between the outflow momentum flux F and protostellar 
bolometric luminosity L DO i as 

log (F/A/ekms-V -1 ) = -5.6 + 0.9 log (L bol /L ). (10) 

In Figure 1141 we plot ted the same correlatio n line (eq. I10|l 
as shown in Fig. 5 of lBontemps et all (|l996h . In the figure, 
Class protostars are distributed in the upper left area to- 
wards the correlation line (solid line), while Class I (and II) 
protostars are roughly distributed along the correlation line. 
Several Class I protostars are distributed in the lower left 
or lower right area towards the line. Because the lifetime 
of Class I protostars is longer than that of Class proto- 
stars, the number of Class I protostars plotted in the figure 
is greater than that of Class protostars. The horizontal dis- 
persion (or luminosity dispersion) is caused by the different 
evolutionary stages of the protostar (or different luminosi- 
ties), while the vertical dispersion (or momentum flux), is 
likely attributed to different cloud parameters. Protostars 
formed in clouds with stronger magnetic fields or slower ro- 
tations are distributed in the lower right area. Thus, cloud 
parameters can be expected from observations of outflow 
momentum and protostellar luminosity. 

In the figure, some Class protostars are distributed 
in the (upper) left area, which indicates that the protostel- 
lar outflow emerges in a very early phase of Class stage, 
or prior to protostar formation. When protostellar outflow 
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appears in the later evolutionary stage, protostars are dis- 
tributed only in the lower right area towards the correlation 
line, as shown by model 6. A comparison of the evolutionary 
tracks in the upper panels with open circles in lower panels 
reveals that no circle is plotted in the region of Lboi <, 0.3 Lq 
and F > 5 x 10 _5 Mq kms -1 yr -1 because such protostars 
have an age of t ps < 1000 yr. If a protostar located in such 
an area is observationally confirmed, it can be identified as a 
very young object (< 1000 yr), through which the very early 
stage of star formation can be investigated. 



4.2 Momentum Flux vs. Envelope Mass 

Figure [15] shows the momentum flux against the envelope 
mass, in w hich the same co r relati on line is plotted as that 
in Fig. 6 of iBontemps et all (|l996l ). 

log {F/Mq kms" 1 ) = -4.15 + 1.1 log (M onv /M ). (11) 

In the figure, the evolutionary track for each model moves 
from near the upper boundary to the lower left. In all mod- 
els except models 6 and 7, outflow appears prior to pro- 
tostar formation. The envelope mass rarely decreases prior 
to protostar formation because the protostar cannot gain 
the mass of the infalling envelope. Even before protostar 
formation, the infalling envelope slightly decreases because 
part of the gas (~ 0.01 Mq) has accreted onto the first core 
(i.e. the rotating disk). In contrast, the outflow momentum 
flux continues to decrease from its emergence (Figure [15)) . 
Therefore, prior to protostar formation, the momentum flux 
moves vertically downward nearly maintaining the initial en- 
velope mass. Then, following protostar formation, the enve- 
lope mass rapidly decreases (Fig. [$]), and the evolutionary 
track moves to the lower left. On the contrary, the evolu- 
tionary tracks for models 6 and 7 start from different points 
(near the center of the figure) in Figure 115b because the 
outflow appears long after the protostar formation during 
which time the envelope mass continues to decrease. 

Models with initially different cloud masses (models 9 
and 10) are also plotted in the left panels in Figure [TBI Be- 
cause the host cloud for these models is more massive than 
that in other models, the evolutionary tracks start from 
M env ~ 1.6 M Q and 2.1 Mq, respectively (Table [3]). How- 
ever, the evolutionary tracks for these models have the same 
trend as that in other models. Thus, different cloud mass 
produces no qualitative differences in evolutionary tracks. 

The outflow momentum fluxes at every thousand years 
plotted in the lower panels of the figures are mainly dis- 
tributed near the correlation line, which indicates that 
the outflow momentum flu xes derived in this cal culation 
agree with the observation of lBontemps et all (| 19961 ). More- 
over, we compared these panel s with the observa tions of 
lHatchell et~all (|2007l . Fig. 4) and lCurtis et al.l (|2010l . Fig. 6), 
and confirmed good agreement between our results and ob- 
servations. In addition, Figure[TS]shows that an outflow with 
F > 3 x 10~ 4 Mq kms -1 yr" 1 indicates a very early phase 
of the outflow (before protostar formation). 



4.3 Evolutionary Relation between Outflow and 
Cloud Core 

With a substantial amount of observational data, 
IBontemps et al.l (|l996T ) expected that the decrease in out- 
flow momentum flux and envelope mass is an evolutionary 
effect independent of protostellar luminosity, or protostel- 
lar mass. They also demonstrated that the outflow mo- 
mentum flux (F) and envelope mass (M anv ) are propor- 
tional to the bolometric luminosity as F cx L^'^ ±0 ' 2 and 
M onv tx LjJ'J 6 2 , respectively. Then, to r emove any lumi- 
nosit y dependence, they created a diagram jBontemps et al.l 
1 19961 Fig. 7), in which the outflow efficiency (F ■ c/Lhoi\ di- 
mensionless, where c is the speed of light) is plotted against 
M e nv/£b 6 i. in their diagram, outflow momentum fluxes for 
Class and I protostars are separately distributed. Class 
protostars are widely distributed, while Class I protostars 
are clustered in a narrow area. They concluded that out- 
flow momentum flux depends on only one basic parameter 
related to M env rather than separately on age and stellar 
mass. 

Our results support their conclusion. To compare obser- 
vations with our results in more detail, we used our calcu- 
lation d ata to create the same diagram (Fig. I16p as that in 
Fig. 7 of lBontemps et alj (|l996l ^l . in which F-c/L hoi (= j/b) is 
plotted against (M e nv /Lboi)° 6 (= xb) for m odels 2, 3, 4, 7, 8 
and 9. As shown in iBontemps et al.l jl996|), the Class and 
I protostars are separately distributed in Figure [TS] Class I 
protostars are clustered in a narrow area of 0.03 < xb < 0.2 
and 10 < t/B < 120 for our results and app eared in an area 
of 0.0 1 < xb < 0.2 and 10 < y B < 600 for IBontemps et al.l 
(|l996h . In addition, Class protostars are widely distributed 
in th e range of xb > . 07 an d 30 < j/b < 10 4 for our results, 
while IBontemps et al.l (|l996l ) showed results of xb > 0.2 and 
200 < j/b < 10 4 . The different distribution between Class 
and I protostars is caused by different conditions of outflow, 
or different evolutionary stages. As shown in Figure [TOl the 
outflow momentum inside the host cloud begins to decrease 
by the end of the Class stage. Thus, during the Class 
I stage, no powerful outflow is driven by the circumstellar 
disk. That is, the protostellar outflow cannot gain additional 
momentum during the Class I stage. However, the outflow 
momentum is nearly conserved and the outflow propagates 
into the interstellar space keeping the momentum acquired 
during the Class stage. It should be noted that the out- 
flow momentum gradually decreases over a lengthy period 
because it interacts with the interstellar medium. Therefore, 
the outflow momentum hardly changes during the Class I 
stage, or the momentum driven phase. In addition, the main 
accretion phase has already ended by the time of the Class 
I stage, as shown in Figure 1111 Thus, the protostellar lu- 
minosity is mainly supplied by the Kelvin-Helmholtz con- 
traction rather than by release of accretion energy. Because 
the Kelvin-Helmholtz timescale is as long as ~ 10 6 yr just af- 
ter the mass accretion ceases, protostellar luminosity change 
is minimal during this phase. As a result, the evolutionary 
track during the Class I stage rarely moves in the vertical 
direction because both F and Lboi are hardly changed. In 
addition, the protostellar luminosity is roughly proportional 
to M onv cx Lboj ( e -g- Lhoi °c M^ 7 ) during the Class I stage, 
as described in the upper panel in Figure [131 Thus, the evo- 
lutionary track also rarely moves in the horizontal direction. 
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As a result, the evolutionary track remains in a small area 
during the Class I stage. 

On the other hand, just after protostar formation, both 
the outflow momentum (Fig. 1101) and luminosity (Fig. I13p 
increase. Next, during the Class stage, protostellar lumi- 
nosity turns to decrease after the outflow momentum turns 
to decrease. The infalling gas first accretes onto the circum- 
stellar disk. Part of the accreted gas in the circumstellar disk 
is blown away by the protostellar outflow, which is powered 
by the accretion from the infalling envelope. Thus, as the 
gas accretion onto the circumstellar disk weakens or the in- 
falling envelope is depleted, the outflow weakens, and its 
momentum inside the host cloud begins to decrease. On the 
contrary, protostellar luminosity is mainly related to gas ac- 
cretion from the circumstellar disk. The accreted gas from 
the infalling envelope remains in the circumstellar disk for a 
short period. Thus, protostellar luminosity does not decrease 
just after depletion of the infalling envelope because proto- 
star luminosity is caused by accretion from the circumstellar 
disk. Therefore, near the end of the Class stage, the outflow 
momentum rarely increases (or the outflow momentum flux 
F begins to decrease), then the protostellar luminosity L de- 
creases some time later. As a result, the evolutionary track 
of F/L 00 \ moves downward in Figure [T6l because the numer- 
ator first decreases. In addition, at this stage, protostellar 
luminosity decreases, or increases, minimally while the in- 
falling envelope rapidly decreases, as shown in Figure IT3T a). 
Thus, the ratio of the envelope mass to protostellar luminos- 
ity Menv/ibol tends to decrease, and the evolutionary track 
moves towards the left. Therefore, during the Class stage, 
the evolutionary track moves to the lower left as shown in 
Figure 1161 It should be noted that the evolutionary track 
does not move widely left because the protostellar luminos- 
ity by the Kelvin-Helmholtz contraction dominates that by 
accretion in the later Class I stage and the protostar is less 
dark. 

In summary, reflecting the rapid attenuation of outflow 
(momentum flux) during the Class stage, the evolutionary 
track moves in a wide range, and the protostar is widely dis- 
tributed, as shown in Figure [16] During the Class I stage, 
the outflow is in a momentum driven, or snow-plough, phase 
without additional momentum. Thus, the evolutionary track 
remains in a small area. Therefore, the different distribution 
of protostars on the diagram between Class and I proto- 
stars is caused by the outflow condition such that outflow is 
powered by mass accretion during the Class stage and is 
in a momentum driven phase during the Class I stage. 



5 DISCUSSION 

5.1 Dependence of Host Cloud Mass 

In this study, we investigated the cloud evolution with dif- 
ferent cloud parameters of the magnetic field Bo, rota- 
tion rate f2o and initial cloud mass M c \. As listed in Ta- 
ble [1] we adopted a wide range of magnetic field strength 
(7.8 x ltrVC < Bo < 7.4 x 10~ 5 fiG) and rotation 
rate (1.0 x 10~ 14 s _1 < fi < 2.1 x 10~ 13 s _1 ), while we 
only adopted three different initial cloud masses of M c i = 
1.05 M , 1.6 Mq and 2.1 Mq. The purpose of this study is 
to associate the properties of the protostellar outflow with 



the envelope mass. Thus, we may have to investigate the 
cloud evolution with a more wide range of initial cloud mass. 
However, with an initially massive cloud core, we cannot 
calculate the cloud evolution until the infalling envelope is 
depleted because it needs a huge computational time. In this 
study, we could calculate the cloud evolution until the mass 
of the infalling envelope decreases to < 10% of the initial 
cloud mass for models with M c \ = 1.05 Mq (models 1-8), 
while 52% for model 10 which has the host cloud mass of 
M c i = 2.1 Mq. Thus, it is considerably difficult to calculate 
the cloud evolution until Class I and II stages with a massive 
host cloud. However, as a result of the calculation, formed 
protostars have a mass of 0.24 — 0.71 Mq, as summarized 
in Table [2] T he initial mas s function seems to have peak at 
£ 1 Mq re.g. lKroupall2001f ). and the cloud mass function in 
various star formi ng regions has a peak around 1 Mq (e.g. 
I Andre et al]|201Ct ). Thus, we believe that we investigated 
the evolution of most typical clouds observed in star form- 
ing regions, or the evolution of typical stars, in this study. 
On the contrary, a massive protostar formed in a massive 
cloud (^> 1 Mq) shows a more massive, or larger momen- 
tum, outflow. Although such massive outflows are small in 
number, they are considered to be preferentially observed. 
In this subsection, we discuss the dependence of the initial 
cloud mass on the protostellar outflow. 

5.1.1 Outflow Energy 

In the left panel in Figure [17] , the outflow energies E ut , 
which is defined as 

E out = i J pout Wout dv, (12) 

where p ou t and twt are the density and velocity of outflow 
at each point, respectively, is plotted against the normalized 
envelope mass for models 3, 4, 7, 8, 9 and 10. The outflow 
energy has a peak around ~ 10 36 — 10 37 erg, and begins to 
decrease during the Class I stage. The outflow energies de- 
rived in o ur calculation are co mparable to the observational 
estimates. ICurtis et all (|2010l ) observed 45 outflows in the 
Perseus molecular clouds and showed that the outflow en- 
ergies lie in the range of 10 35 erg < E ou t < 10 37 erg except 
for a highly powerful anomaly. On the contrary, a massive 
proto star seems to drive a more powe rful outflow (|Wu et al. I 
12004? ) . Recently. iMotoei et ail (|201lh showed that the out- 
flow around a young massive protostar that is embedded in 
~ 200 Mq envelope has the kinetic energy of > 10 46 erg. 

The left panel in Figure [17] indicates that the outflow 
kinetic energy strongly depends on the initial cloud mass, 
while it weakly depends on the cloud parameters of the mag- 
netic field and rotation rate. The outflow kinematic energy 
for model 10 is about ten times larger than that for model 
4 that has the same cloud parameters of magnetic field (6) 
and rotation (w) as in model 10 but the different host cloud 
mass (Mci = 2.1 M for model 10, and 1.1 Mq for model 4). 
The outflow kinetic energy (and momentum) continues to 
increases until the infalling envelope almost halves (e.g. dur- 
ing the Class I stage: Fig. I10[) . A massive cloud has a longer 
duration of the Class stage during which a sufficient mas- 
sive infalling envelope can give power to drive the outflow. 
Note that the outflow energy largely dissipates inside the 
host cloud by the interaction between outflow and infalling 



© 0000 RAS, MNRAS 000, 000-000 



Evolution of Protostellar Outflow 13 



envelope, in which the outflow loses its energy radiatively 
at the shock. Note also that since the mass of the infalling 
envelope depends on the initial cloud mass, outflows are ex- 
pected to have different energies among models with differ- 
ent cloud mass. As a result, the protostellar system formed 
in a more massive envelope can drive the protostellar out- 
flow for a longer duration, and outflow in such a system can 
acquire a larger kinematic energy, or momentum. Thus, it is 
expected that outflow with a larger kinematic energy tends 
to be observed in a massive infalling envelope. However, it is 
also expected that such powerful outflows are not ubiquitous 
due to the rarity of a massive cloud core (M c i 3> 1M©). 

5.1.2 Outflow Kinematic Luminosity 

In the right panel in Figure 1171 the outflow kinematic lumi- 
nosity L m , which is defined as 

i m = f^, (13) 

tout 

is plotted against the protostellar luminosity for every 
3000 yr after the protostar formation. The evolutionary 
track for model 4 (the blue solid line) is also plotted. 
In the figure, protostars are distributed in the range of 
0.5 L© < Lboi < 501/0 and 0.005 L© < L m < 0.07 L . 
In the same range of the bolometric luminosity, almost the 
same distribution of protostars is seen i n the observations 
jCabrit fc Bertoutlll992l ; IWu et al. 1120041 ). 

Both our results and observations indicate that the me- 
chanical luminosity weakly depends on the protostellar bolo- 
metric luminosity. In the panel, the evolutionary track of 
L m is much shallower than L m oc Lboi; it moves horizontally 
during the main accretion phase and moves to the lower left 
after the main accretion phase. This means that the out- 
flow driving is not directly related to protostellar properties 
such as luminosity, mass and age. The protostellar proper- 
ties are determined by the accretion process onto the proto- 
star from the circumstellar disk. On the other hand, outflow 
properties are determined by the accretion process onto the 
circumstellar disk from the infalling envelope. Thus, the pro- 
tostellar outflow is related to the envelope mass, while the 
protostellar luminosity is related t o the circumstel l ar dis k 
mass. Therefore, as pointed out by iBontemps et "all (|l996T ). 
it seems that the protostellar outflow is not primary related 
to the protostellar luminosity. This implies that the proto- 
stellar outflow is not powered by the radiation of the pro- 
tostar. In addition, there is no significant difference in the 
peak value of the mechanical luminosity among models with 
different host cloud masses. This indicates that the outflow 
is powered by the accretion at a constant rate because the 
mechanical luminosity, which is the outflow energy divided 
by its lifetime, is almost constant during the main accretion 
phase. 

In summary, the acquisition rate of the outflow energy 
or momentum is independent of the initial cloud mass, while 
it slightly depends on the cloud parameters of the magnetic 
field and rotation rate. The protostellar outflow gains its 
energy or momentum at a constant rate during the main ac- 
cretion phase. Since a massive cloud has a longer period of 
the main accretion phase that almost corresponds to dura- 
tion of the Class stage, the outflow appeared in a massive 
cloud core can possess a larger energy and momentum. In 



other words, a massive star shows a more powerful outflow 
because a relatively massive star forms in a massive cloud. 
Therefore, the outflow (peak) energy, or momentum, reflects 
its host cloud mass. 



5.2 Cloud Parameters 

In §4] we have shown that the outflow in the Class stage 
is e ssentially different fro m that in the Class I stage (see 
also lBontemps et al.lll996T l. The circumstellar disk forcefully 
drives outflow during the Class stage, while the driving 
force from the circumstellar disk ceases and outflow is in 
the momentum-driven (slow-plough) phase during the Class 
I stage. As a result, Class and I protostars are distributed 
in different regions in Figure 1141 in which the outflow mo- 
mentum flux is plotted against the protostellar bolometric 
luminosity. 

Observations often show a large scatter of the outflow 
momentum flux even am ong protostars whic h are in the 
same evolutionary stage |Hatchell et alj I2007T 1 . Such scat- 
ter may be attributed to variations of the cloud parameters 
such as the strength of the magnetic field and rotation rate. 
In Figure 1141 the protostar formed in a cloud with strong 
magnetic field (e.g. model 5) or slower rotation (model 6 
and 7) has a relatively smaller momentum flux. The outflow 
momentum in these models is about one order of magnitude 
smaller than the other models. A scatter of the mechani- 
cal luminosity is also seen in the right panel in Figure 1171 
Thus, it is natural that the observation shows a scatter of 
outflow pro perties because the host cloud has different cloud 
parameters (|Casellj|2063 : ICrutcherlll999l ). 

5.3 Protostellar Mass and Star Formation 
Efficiency 

The protostellar outflow reduces the star formation effi- 
ciency because it ejects materials from the host cloud into 
the interstella r space. The magnetic field and rotation drive 
the o utflow jBlandford fc Payne! Il982l ; lUchida fc Shibatal 
Il985l ). Thus, it is expected that both magnetic field and 
cloud rotation control the outflow and star formation effi- 
ciency. As summarized in Tabled even in clouds having the 
same mass, different magnetic field strengths and rotation 
rates bring different results. For example, a half of the cloud 
mass is ejected by the protostellar outflow for model 2, while 
only 10% of the host cloud mass is ejected for models 5 and 
6, in which the magnetic field strength for models 5 and 
6 is stronger than that for model 2. As described in t)3.3l 
both a moderate strength of the magnetic field and rapid 
rotation rate are necessary to drive a powerful outflow. Nei- 
ther weak magnetic field nor slow rotation can drive a pow- 
erful outflow. In addition, in a strongly magnetized cloud, 
the magnetic braking effectively transfers the disk angular 
momentum and delays the formatio n of a sufficiently larg e 
disk that is the driver of the outflow jMachida et aLlfeoilbl ) . 
Thus, too strong magnetic field also suppresses a powerful 
outflow driving. 

In models with the same cloud mass of M c \ — 1.05 Af© 
(models 1-8), the protostars have their mass in the range 
of M ps = 0.24 — 0.71 Af©, which corresponds to the star 
formation efficiency of e = 0.23 — 0.68. As seen in Table [2j 
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the star formation efficiency is high in a cloud with stronger 
magnetic field or slower rotation. In such a cloud, a relatively 
small circumstellar disk drives a relatively weak outflow that 
cannot eject a large fraction of cloud mass from the host 
cloud, and a relatively massive protostar finally forms. Thus, 
the outflow efficiency is inversely correlated with the star 
formation efficiency. 

The outflow efficiency is also inversely correlated with 
the protostellar luminosity in the main accretion phase. As 
described in Table[2] at the end of the Class stage, the more 
luminous protostar has the less massive outflow (models 5, 
6 and 7) and vice versa (models 1,2 and 3). Figure Q3] also 
shows that the luminous protostars have smaller outflow mo- 
menta. The protostellar luminosity is related to the outflow 
efficiency through the efficiency of the angular momentum 
transfer in the circumstellar disk. In a strongly magnetized 
cloud, the disk angular momentum is effectively transferred 
not by the outflow but by the magnetic braking that pro- 
motes the mass accretion from the circumstellar disk onto 
the protostar without mass ejection, and a relatively mas- 
sive protostar forms. Since the accretion luminosity is pro- 
portional to both the accretion rate and protostellar mass, 
both the higher mass accretion and massive central object 
makes the protostar luminous. 

On the contrary, with a relatively weak magnetic field, 
the angular momentum transfer is inefficient with the weak 
magnetic braking. In such clouds, the strong disk rotation 
drives the powerful outflow suppressing the rapid (or effec- 
tive) mass accretion onto the protostar. Therefore, the less 
luminous protostar tends to have a more powerful outflow. 

Note that this trend may not be applicable among 
clouds with initially different masses because a massive pro- 
tostar can form with a powerful outflow when the initial 
cloud has a sufficient mass, irrespective of the magnetic field 
strength and rotation rate. However, we can observe these 
trends in some star forming regions where low-mass stars 
born in clouds with similar mass. We need a more detailed 
observation of the core mass function in various star form- 
ing regions to confirm the relation between the protostellar 
luminosity and outflow properties. 

5.4 Molecular Outflow Model 

A very young protostar sometimes shows two, or more, 
different types of flow: high- velocity (or optical) jet and 
low-velocity (or molecular) outflow. The mole cular out- 
flow i s frequently observ e d around a protostar (|Wu et al. I 
|2004 IZhang et al] 120051 ; lHatchell et all l2007h . while the 
optical (or the high-velocity mole cule) jet is occasion- 
ally observed in molecular outflow jMundt fc Fried! 1 19831 ; 
Richer et all 1 19921; iMtchell et all Il994l ; lArce fe Goodmanl 



2002; iVelusamv et al.1 120071 ). Observations indicate that a 



high-velocity jet is enclosed by the low-velocity outflow. 
Since molecular outflow has a large amount of mass (e.g. 
iDownes fc Cabritl 12003] ; IStoiimirovic et al.1 120061 ). it affects 
the star formation effic iency and greatly c ontributes to the 
star formation process l|Nakano et alj[l995l ). However, in ob- 
servations, it is difficult to specify the driving mechanism 
for molecular outflow because we cannot directly observe 
the driving region that is embedded in a dense cloud core. 
Thus, many authors have tried to theoretically clarify it. 
Historically, the entrainment mechanism has been pro- 



posed to explain the molecular outflow driving, in which 
one imagined the low- velocity outflow entrained by a high- 
velocity jet. In other words, the primary jet injects its mo- 
mentum into the surrounding gas resulting in molecular 
outflow. This mechanism was simply formulated with ob- 
servations of high- velocity well-colli mated jets embedded i n 
the wide-angle low- velocity outflow (|Mundt fc Friedlll983l V 
There are a vari ety of entrainment models such as th e 
wind-driven shell (IShu et aI.lll99ll;lMatzner fc McKedl 19991) . 
jet-dri ven bow shock ( Raga fc Cabritj ' 1993al; Raga et al] 



1993bl). and jet-driven t urbulent ( Canto fc Ragal 1991 



Lizano fc Giovanardilll995r ) models. It seems that some mod- 



els are possible to explain a number of observations of molec- 
ular outflow and some is difficult to explain general feature 
of ou tflows (|Cabrit et al] 1 19971 ; [Lee et al][2000l ; lArce et al~l 
l2007t ). Around the protostar, the molecular outflows are 
frequently observed, while the high-velocity component is 
rarely observed. Thus, we have to posit the undetected jet 
to the model molecular outflow with the entrainment mech- 
anism, although the high-velocity jet may be invisible by 
chance. In addition, to explaining observations, one can ad- 
just jet and outflow models by changing a number of arti- 
ficial, or ambiguous, parameters such as the jet speed and 
momentum conversion efficiency between jet and outflow, 
environment density, strength of magnetic field, etc. 

Aside from the entrainment model, there is an entirely- 
different appr oach for theoreti cally investigating the molec- 
ular outflow. iTomisakal (|2002[ ) calculated the evolution of 
the whole molecular cloud without any artificial setting and 
reproduc ed the low-velocity ou tflow in the collapsing cloud 
(see also ITomisakal 1 19981 . |2000| V In his study, the first core 
forms before the protostar formation and directly drives the 
low-velocity outflow. He also reproduced the high-velocity 
jet which is driven near the protostar. Then, both low- 
and high-velocity flows appeared in the c ollapsing cloud 



have been confirmed by many a uthor s (e.g. iMachida et al 



20041; IMatsumoto fc Tomisakal [2004; 



Bancrice & Pudritz 



2006; iHennebelle fc Fromang 2008; Commc rcon et al 



t al.ll2010al ; lBurzle et al]|201ll ; ISeifried et al 
20121 ) 



Price et al 



2010 



2012; 



The high-velocity jet and low-velocity 
outflow are naturally reproduced without any artificial set- 
ting in cloud collapse calculations, while the jet is artificially 
input to entrain the molecular outflow in entrainment mod- 
els. However, in the cloud collapse calculations, the authors 
could not investigate a long-term evolution of jet and outflow 
because they had to resolve the protostar i t self that is the 
driver of the jet jBaneriee fc Pudritzl 120061 ; IMachida et al] 
l2008bl ; iTomida et alj|2012h . 

In this study, at the expense of spatial resolution around 
the protostar, we could investigate the evolution of molec- 
ular outflow and protostellar system for ~ 10 yr with sink 
treatment. However, no high-velocity jet appears because 
we did not resolve the region in the proximity of the pro- 
tostar (r < 1 AU) w here the high- velocity jet is driven 
|Machida et alj|2008bh . Thus, we could not investigate the 
effect of the jet on the molecular outflow. On the con- 
trary, since we resolve t he magnetically i nacti v e region inside 
the c ircumstellar disk (|Machida et al] 120071 ; ITomida et al] 
l2012h , we could precisely investigate the low- velocity outflow 
driven by the circumstellar disk. Note that the low-velocity 
outfl ow is launched outside of the magne tically inactive re- 
gion jlnutsuka et al]|2010l ; Ilnutsukall2012l ). 
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In [|3]and[4] we have shown that our results agree well 
with observations of molecular outflows; the outflow momen- 
tum (flux) and energy are comparable to the observations. 
Since we have not focused on the effect of a high-velocity 
jet, we do not exclude the entrainment mechanism. In re- 
ality, a part of the infalling gas might be entrained by the 
high-velocity jet. However, our results suggest that the ob- 
served molecular outflows can be explained only by the flow 
directly driven by the circumstellar disk. In contrast to the 
entrainment models, the direct cloud collapse simulations 
explain the molecular outflow well without any fine-tuning of 
model parameters. Although further simulations with high 
spatial resolutions would be useful to examine the effect of 
the high-speed jet, our current simulations mostly explain 
the available observations of the low-velocity outflow. 



down to ~ 50%, whose exact value depends on the cloud 
parameters such as the magnetic field strength and rota- 
tion rate. A stronger magnetic field excessively transfers the 
angular momentum and forms a relatively small circumstel- 
lar disk that drives a relatively weak outflow. A weak mag- 
netic field and slow cloud rotation also weaken the outflow 
driving force. On the contrary, the clouds with moderate 
magnetic fields and rotation, bo = 0.05 — 0.4 (fj, ~ 3 — 10) 
and /3o = 0.01 — 0.04, show considerably powerful outflows. 
This para meter range agrees with th e observational esti- 
mates (e.g. lCrutcherll 19991 ; [Caseilill2002l ). Therefore, the pro- 
tostellar outflow should determine the final stellar mass and 
significantly affect the early evolution of the low-mass pro- 
tostars. 



6 SUMMARY 

In this study, we calculated the cloud evolution from the 
pre-stellar core stage until almost all the envelope gas dis- 
sipates by the protostellar outflow. In the collapsing cloud, 
the first core forms prior to the protostar formation and 
evolves into the circumstellar disk. The low- velocity outflow 
is first driven by the first core, then driven from the outer 
part of the disk which is magnetically active. The outflow 
driving region extends as the circumstellar disk grows. Be- 
fore the outflow breaks out of the cloud, the outflow propa- 
gates along hourglass-shaped magnetic field lines that open 
up with increasing the distance from the equator. As a re- 
sult, the outflow has a wide-opening angle in a cloud scale. 
After the outflow penetrates the cloud, it propagates along 
the interstellar magnetic field lines. The head of the outflow 
travels over ~ 10 s AU in 10 5 yr. In contrast, the horizon- 
tal extension of the outflow is limited by the cloud scale. 
The width of the outflow reflects its host cloud size. The 
outflow thus extends only in the vertical direction, and its 
collimation gets improved in this stage. 

The properties of the calculated outflows such as the 
outflow momentum, energy and mass agree well with those 
of observed outflows. Our simulations also explain the phys- 
ical structure of the observed outflows. These support the 
picture that the low-velocity molecular outflows which are 
frequently observed around protostars are directly driven in 
the circumstellar disk, or the first core. The entrained flow is 
not necessary to explain observations, though some amount 
of gas might be entrained by the high- velocity component. 

Our calculations show the same correlations between 
the outflow momentum flux, protostellar luminosity and en- 
velope mass as in observations. These correlations differ be- 
tween Class and I protostars, which is explained with the 
different evolutionary stages of low-velocity outflow. In the 
Class stage, the sufficient gas accretes onto the circumstel- 
lar disk, and the outflow powered by the accretion is driven 
from the disk. The outflow momentum is continuously sup- 
plied from the disk, or the accreting gas, during this stage. 
However, the outflow gradually ceases as the infalling enve- 
lope gets depleted. In the Class I stage, the outflow hardly 
acquires its momentum from the accreting gas and enters 
into the momentum-driven snow-plough phase. 

The protostellar outflow ejects half of the cloud mass 
from the host cloud and limits the star formation efficiency 
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Table 1. Model parameters 
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Table 2. Calculation Results 
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8 


0.25 


0.12 


0.31 


12.5 


3.8 xlO 4 


0.33 


0.12 


0.5 


2.5 


l.lxlO 5 


0.34 


0.13 


0.48 


0.09 


9 


0.23 


0.30 


0.46 


4.1 


8.9 xlO 4 












0.27 


0.37 


0.53 


0.49 


10 






















0.49 


0.28 


0.30 


1.10 



* E.O.C means the end of calculation. 



Table 3. Protostar formation and Outflow Emergence epochs 



Model t ps ,o [yr] i ou t,o [yr] 



1 


7.5714 


X 


10 4 


7.2963 


X 


10 4 


2 


7.4285 


X 


10 4 


7.0360 


X 


10 4 


3 


7.3127 


X 


10 4 


7.0466 


X 


10 4 


1 


7.5860 


X 


10 4 


7.1786 


X 


10 4 


5 


8.2225 


X 


10 4 


7.8861 


X 


10 4 


6 


6.8280 


X 


10 4 


8.1441 


X 


10 4 


7 


7.0053 


X 


10 4 


7.3155 


X 


10 4 


8 


9.2215 


X 


10 4 


7.5304 


X 


10 4 


9 


9.3039 


X 


10 4 


7.5246 


X 


10 4 


10 


8.8585 


X 


10 4 


8.0162 


X 


10 4 
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Figure 1. Structure of protostellar outflow at f ps = 1.447 X 
10 5 yr after the protostar formation with different scales for model 
3. The box scale is described in each panel. The elongated blue 
and red structure in panels (a)-(d) mean the protostellar outflow 
inside which the gas is outflowing from the central region. The 
central sphere with blue surface in panels (a) and (b) corresponds 
to the host cloud. The colour on the cutting plane inside the host 
cloud means the density distribution. The pscudo disk, infalling 
envelope and rotating disk are indicated by an arrow in each 
panel. The thick arrows with black and white colour in panel (e) 
are velocity vectors in the outflowing region. The magnetic field 
lines inside the z > outflowing region are plotted by streamlines 
in panels (d) and (e). 
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Figure 2. Time sequence of the protostellar outflow in a cloud 
scale for model 3. The protostellar outflow is represented by yellow 
and orange iso- velocity surface inside which the gas has a positive 
radial velocity (v r > 0). The inner host cloud region of 0° < 6 < 
90° and 0° < <j> < 90° is plotted. The colour on each wall is the 
density distribution on the cutting plane. The elapsed t and t ps 
are described in each panel. 
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Figure 3. Evolution of the protostellar outflow in a large scale 
for model 3. fn each panel, the central grey sphere is the host 
cloud with a radius of r = 6.1 X 10 3 AU. The scale of each panel 
is different. The outflow colour means the 2-component of the 
outflow velocity, v z . 
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-6000 6000 



x [AU] 

Figure 4. Shape of the outflow in the region of z > with 
different epochs for model 3. Each epoch (the elapsed time after 
the protostar formation) is described in the upper right region. 
The angle between two arrows in the left panel mean the outflow 
opening angle at each epoch, while that in the lower right panel 
is the opening angle in a cloud scale. The grey arrows in the lower 
right panel is the evolutionary direction of the outflow opening 
angle in a cloud scale. The red circle corresponds to the host 
cloud. 
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Figure 5. Evolution of the density (each upper panel) and out- 
flow momentum (each lower panel) distribution on the y = 
cutting plane in the host cloud scale for model 3. The white line 
in each panel corresponds to the boundary between the outflow- 
ing gas and infalling envelope. The host cloud (r = R c i) is plotted 
by the red circle. Arrows are velocity vectors (only the velocity 
vectors in the outflowing region are plotted in each lower panel). 
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Figure 6. Same as in Fig. [5] but in the rotating disk scale. 
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Figure 7. (a) Mass of protostar, disk, outflow and envelope, (c) 
total and protostellar luminosities (left axis) and mass accretion 
rate (right axis) and (e) outflow momentum (left axis) and energy 
(right axis) for model 3 are plotted against the elapsed time after 
the cloud begins to collapse. The outflow mass in the host cloud 
(r < R c \) and that in the whole computational domain are plotted 
in panel (a). The momentum and energy of outflow in the host 
cloud (r < i? c i) and those in the whole computational domain 
are plotted in panel (e). The same physical quantities are plotted 
against the envelope mass normalized by the initial host cloud 
mass in the right panels (b), (d) and (f). 
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Figure 8. Time evolution of the envelope mass for models 1- 
8. The left axis is the envelope mass normalized by the initial 
cloud mass, and right axis is the cloud mass in unit of Mq. The 
evolutionary stages (Class 0, I and II) are described according to 
the classification in i|3.1l The elapsed time t = 5tg q and lOtgf rj 
are plotted by dash lines, where tg q is the free fall timescale of 
the initial cloud. The relation M cnv oc t~ 2,5 is also plotted by the 
thick line. 
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Figure 9. Mass of protostar, disk and outflow for models 2, 4, 
5, 6, 7 and 8 are plotted against the envelope mass normalized by 
the initial cloud mass, respectively. In each panel, the blue solid 
lines are the outflow mass in the whole computational domain, 
while the blue broken lines are the outflow mass inside the host 
cloud r < R c \. 



© 0000 RAS, MNRAS 000, 000-000 



28 M. N. Machida, & T. Hosokawa 



Class | Class f ! Class II 

T 1 i i 1 1 1 1 1 < 1 1 1 1 1 1 — 




Figure 10. Outflow momentum for models with different initial 
magnetic field strengths (upper panels) and different initial cloud 
rotation rates (lower panels) against the envelope mass normal- 
ized by the initial cloud mass (left panels) and the time after 
the protostellar outflow emerges (right panels). In each model, 
the outflow momentum in the whole region (solid line) and that 
inside the host cloud (broken line) are plotted. 
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Figure 11. Protostellar luminosity (blue line; left axis) and mass 
accretion rate (black and red lines; right axis) for models 2, 4, 5, 
6, 7 and 8 are plotted against the elapsed time after the protostar 
formation. The red line in each panel is the mass accretion rate 
averaged every 1000 yr. 
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Figure 12. Density (colours and contours) and velocity (arrows) 
distributions around the protostar on the equatorial plane for 
model 8. The elapsed time is described in the upper left corner. 
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Figure 13. Protostellar luminosities for models 1, 2, 3, 5, 6 
and 7 are plotted against the normalized envelope mass (upper 
panel), protostellar mass (middle panel) and elapsed time after 
the protostar formation (lower panel), respectively. The relation 
of L oc Af^^ 7 (M onv oc L 6 ) is plotted in upper panel. 
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Figure 14. Outflow momentum fluxes (F) for models with dif- 
ferent magnetic field strengths (left panels) and different initial 
cloud rotation rates against the protostcllar luminosity are plot- 
ted in the upper panel. The outflow momentum flux every 1000 yr 
are plotted in the lower panels. The open and filled circle mean 
the momentum flux during the Class and I stages, respectively. 
The plot ted ranges and s olid l ine in each panel is the same as 
Fig. 5 of lBontemps et al.l lll996l 1. 
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Figure 15. Same as in Fig. I14l but against the envelope mass. 
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Figure 16. Momentum flux divided by the protostellar luminos- 
ity against the (M env / M )/(L bol /L sun ) - 6 . 
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Figure 17. Outflow energy -Eout against the normalized enve- 
lope mass (left) and outflow kinematic luminosity L m against the 
protostellar bolometric luminosity Z/ DO i ( rl ght) for models 3, 4, 7, 
8, 9 and 10. The relation L m = Lbol is plotted in the right panel. 
The evolutionary track for model 4 is also plotted by the blue 
solid line in the right panel. 
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